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A cooperative  model  has  been  developed  to  describe  ferroelectric  polarization 
in  polyvinylidene  fluoride  (PVDF).  The  molecular  dipoles  within  the  crystal 
are  assumed  to  have  two  or  more  orientations  available  to  them  and  the  lattice 
eiergy  of  a given  orientational  site  is  assumed  to  be  proportional  to  the 
fraction  of  molecules  having  that  orientation.  An  analytical  solution  is 
given  for  the  2-site  model  which  predicts  polarization  hysteresis  tvpical  of 
ferroelectrics . However,  a more  complex  6-site  model,  which  can  be  analysed 
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numerically,  is  needed  to  account  for  observed  infra  red  hysteresis  data 
and  electric-field-induced  x-ray  structural  changes  which  have  been 
reported  for  PVOF.  Although  the  model  is  simple,  rather  complex  behavior 
is  observed  including  a gradual  increase  or  decrease  in  the  remnant 
polarization  with  number  of  cycles  of  electric  field  application.  Though 
the  agreements  with  various  experimental  data  are  good  an  obvious  need  to 
include  kinetic  effects  in  the  model  is  indicated. 
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FERROELECTRIC  POLARIZATION  IN  POLYMERS 

M.  G.  Broadhurst  and  G.  T.  Davis 
National  Bureau  of  Standards 
Washington,  D.C.  20234 

INTRODUCTION 

Previous  work  has  shown  that  a model  of  oriented  dipoles  accounts  quite  well 
for  observed  piezoelectric  and  pyroelectric  behavior  of  amorphous1  and  semicrystal- 
line*1  polymers  The  process  whereby  the  dipoles  become  oriented  is  the  subject  of 
this  pacer.  Here  we  summarize  a model  and  calculated  results  for  ferroelectric 
reorientation  of  molecular  dipoles  in  the  all-trans  polar  crystal  phase  (9  phase 
or  Form  I)  of  polyvinyl idene  fluoride.  By  ferroelectric  we  mean  a material  with 
polar  crystals  where  the  polarization  direction  can  be  reversed  with  an  applied 
electric  field. 

In  the  model,  the  orientations  of  molecular  dipoles  depend  on  the  electric 
field  and  the  history  of  the  field  application.  The  quantities  'Cose>,  <cos~e>, 
and  f(j)  can  be  calculated  where  9 is  the  angle  between  a dipole  and  the  aoolied 
field,  f(e)  is  the  fraction  of  dipoles  at  angle  9 and  the  brackets  <>  indicate  a 
spacial  average.  The  experimental  physical  quantities  of  interest  for  comparison 
with  the  model  are  polarization  (a  measure  of  <cose>),  infra  red  transmission 
intensity3’4  (a  measure  of  <cos“0>)and  x-ray  pole  figure  data5  (a  measure  of  f(9l). 
Since  only  one  adjustable  parameter  is  used  in  the  model,  these  various  comparisons 
permit  an  unusually  thorough  evaluation  of  the  strengths  and  weaknesses  of  the 
model . 

THE  MODEL 

Polyvinyl idene  fluoride  (PVDF)  is  a linear  molecule  with  a carbon  backbone. 
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In  its  crystalline  forms  the  molecule  has  a large  moment  normal  to  the  chain. 
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In  the  6 crystal  phase  the  unit  cell  is  polar  and  the  crystals  form  lamellae  typi- 
cally 10  nm  thick  with  the  molecular  segments  nearly  normal  to  the  large  crystal 
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surfaces.  The  molecular  segments  which  traverse  the  crystal  lamellae  are  inter- 
connected through  the  intercrystalline  amorphous  regions  (See  Fig.  1)  either  by 
tight  folds  at  the  lamellar  surface  or  by  longer  irregular  molecular  segments.  As 
demonstrated  in  the  data  referenced  later  in  this  paper,  the  segments  in  the  crystal 
can  rotate  as  rigid  rods  or  by  a propagating  twist  mechanism  about  their  chain  axis, 
giving  a reorientation  of  the  crystal  moment.  The  shape  of  the  curve  relating 
potential  energy  to  angle  for  individual  chain  rotation  is  not  known,  but  following 
Kepler  and  Anderson’s  observation  that  the  near-hexogonal  unit  cell  cross  section 
will  permit  a six-fold  degeneracy  in  the  orientation  of  the  crystal,5  we  assume 
equivalent  potential  energy  minima,  60  degrees  apart  as  shown  in  Fig.  2.  The 
potential  energy  of  an  orientation  e ^ relative  to  a base  line  is  U-j  and  the 
fraction  of  segments  having  that  orientation  is  f ^ . We  use  the  first  order, 
mean  field  cooperativi ty  assumption6 ,7  that  Ui  =-U0fi  where  U0  is  the  energy 
difference  between  a filled  and  empty  site.  This  assumption  describes  the 
desired  feature  that  if  most  segments  have  a given  orientation  then  that 
orientation  is  the  preferred  one  (having  the  lowest  energy).  The  lattice  energy 
of  the  crystal  is  given  by  sU1-fi  = -U0Ef.j2,  the  dipole-field  interaction  energy 
is  given  by  -m0Elfi  cos®i  (where  the  field  is  assumed  to  be  at  0 degrees),  and  the 
entropy  is  given  by  -kEf^  in  fj . We  do  not  include  dipole-dipole  interactions. 

The  Helmholtz  free  energy  is 

(Ai  - Aq)  = -W-  m0EEfi  cose  i +kT  (1) 

We  can  minimize  the  free  energy  including  the  constraint  that  Ef^  = 1 by  using 
a Lagrange  multiplier  technique.  That  is,  we  solve  the  equations. 


e(Ai  - Aq)/  3fj  -A  3(Ef1)/3fj  = 0. 


(2) 


The  result  is  a set  of  6 equations 


- 2Uq  fi  - m0E  cosei  ♦ kT(l  +in  fi)  - \ » 0 (3) 

We  can  solve  the  simple  2 site  case  analytically8  but  this  6 site  example  must 
be  solved  numerically  as  detailed  elsewhere.8 

RESULTS 

The  result  of  the  calculation  is  that  at  E = 0 for  U0/kT * 2.0118  (an  unusual 
type  of  Curie  point  for  this  model),  the  6-fold-degenerate  lowest  energy  solution 
is  for  one  site  to  be  heavily  populated  and  the  remaining  five  sites  tc  have 
lesser  (but  equal  to  each  other)  populations.  At  Ie)  > 0 , the  degeneracy  is  re- 
moved and  depending  on  E and  the  angle  between  E and  a preferred  site,  a 
particular  distribution  may  become  unstable  and  a new  distribution  form  resulting 
in  a change  in  the  magnitude  and  direction  of  the  crystal  moment.  The  assumed 
mode  of  redi stri bution  of  orientations  is  somewhat  arbitrary  but  the  one  giving 
results  closest  to  experiment  assumes  that  once  a site  i is  unstable,  the  prob- 
ability that,  of  the  remaining  stable  sites,  site  j will  be  favored  is  proDor- 
tional  to  exp(mQE  cos^j/kT).  This  assumption  is  not  to  be  confused  with  the 
equilibrium  assumption  of  Kepler  and  Anderson-  that  the  population  of  site  j 
will  be  proportional  to  this  term.  In  our  case  the  orientations  are  not  at 
equilibrium  and  are  highly  dependent  on  the  initial  distribution  of  orientations 
and  history  of  electric  field  application. 

The  numerical  calculation  is  done  by  varying  mQ E/k T in  equation  3 (after 
dividing  by  kT)  by  small  increments,  recomputing  f ■ and  if  a given  solution  is 
unstable,  redistributing  the  f.  as  described  above.  We  assumed  an  initial 
distribution  of  crystals  such  that  the  crystal  moments  and  applied  field  are 
in  the  plane  of  the  crystal  lamellae  and  the  moments  are  equally  distributed 
every  10  degrees  from  0 to  350  degrees.  From  the  record  of  f^  versus  E we 
calculate  <cose>  , <cos:e>  and  the  fraction  of  reflecting  crystal  planes  at 
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angle  ^ as  a function  of  field  and  history.  To  compare  mQE/kT  to  the  applied  field 
we  assume  a molecular  segment  in  the  crystal  is  10  nm  long.  Such  a segment  contains 
40  repeat  units  of  6.9  \ 10*-'°  Cm  vacuum  moment  each.  The  reaction  field  (in  the 
SDherical  approximation)  enhances  this  moment  by  (ec  + 2)/3  * 5/3  giving  a total 
segment  moment  of  4.6  x 10':3Cm.  Thus,  a 1 MV/cm  field  at  room  temperature  qives 
m^E/kT  3 11.  This  value  is  well  above  the  range  of  linear  dielectric  reSDonse 
typically  encountered  with  molecular  dipoles.  We  used  a lattice  energy  o*  U0/kT  * 3 
for  our  calculations  which  gives  a critical  field  (for  which  the  first  crystals 
switch)  of  about  1 MV/cm.  Remembering  the  UQ  is  the  energy  difference  between  an 
empty  and  filled  site,  we  see  that  the  dipole-field  energy  can  be  much  greater 
than  the  crystal  energies.  Reducing  U0/kT  to  2.5  reduces  the  critical  field  by 
half  but  as  a function  of  reduced  field,  E/E  (critical),  the  calculations  are 
insensitive  to  U0/kT.  For  simplicity,  we  have  ignored  kinetic  effects  in  the 
calcultaions  even  though  they  are  inherent  in  the  model.  For  example,  thermodynamic 
fluctuations  will  allow  a crystal  to  switch  to  a new  oreferred  orientation  at  fields 
below  the  critical  value  and  all  redistributions  of  dipoles  involve  rotations  which 
encounter  local  energy  barriers  - a orocess  which  takes  time. 

COMPARISON  WITH  EXPERIMENT 

The  calculated  polarization  as  a function  of  electric  field  through  1\  cycles 
is  shown  in  Figure  3.  Experimental  curves  obtained  by  cycling  the  field  at  a 
constant  rate  of  change10  are  shown  in  Figure  4.  Both  the  observed  critical 
switching  field  and  degree  of  polarization  (as  a fraction  of  the  total  possible 
polarization)  are  satisfactorily  mimicked  by  the  calculated  curve.  The  quantity 
mgE/kT  * 11  corresponds  to  E » 1 Mv/cm  for  10  nm  thick  crystals  and  for  these 
mixed  phase  samples  (30T<  a crystals,  30*  $ crystals  and  401-  amorphous  ohasel  the 
maximum  polarization  will  be  about  10  uC/cm ‘ . To  account  for  the  reversible 
permittivity  of  the  amorphous  phase  we  have  added  a relative  permittivity  of  15 
to  the  model  which  accounts  for  the  slope  in  curves  of  Fig.  3 upon  decreasing  the 
field. 
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Figure  5 is  a comparison  between  calculated  <os:6>  for  the  model  and  the 
measured  hysteresis  of  the  512  cm'1  IR  absorption  in  PVDF3.  This  absorption  is 
due  tod  phase  crystal  vibrations  polarized  parallel  to  the  CF , dipole  moment. 
Random  orientation  of  single  axes  rotators  with  the  axes  in  the  plane  of  the 
film  corresponds  to  <os-e>  a 0.5  and  23X  transmission.  Complete  alignment 
corresponds  to  <-cos*-e>  = 1,0  and  40i  transmission  and  these  values  are  used 

to  match  the  magnitudes  of  the  ordinates.  Thus  the  magnitudes  and  general  positions 
of  the  two  curves  are  very  similar. 

In  Figure  6 we  show  calculated  and  observed  x-ray  pole  figure  results 
obtained  at  zero  field  after  the  material  has  been  subjected  to  a polarization 
treatment.  The  abscissa  refers  to  the  angle  between  the  bisector  of  the  incident 
and  diffracted  x-ray  beams  and  the  normal  to  the  PVDF  film  (or  equivalently  the 
direction  of  the  poling  field).  The  results  from  the  present  calculation  (dashed 
curve)  predict  the  correct  magnitudes  at  zero  and  large  angles  and  are  less 
satisfactory  than  the  equilibrium  calculation'  (solid  line)  at  small  angles. 

While  the  above  results  support  a conclusion  that  the  essential  features  of 
tne  high  field  Dolarization  in  PVDF  are  described  by  a simple  cooperative  ferro- 
electric model,  the  curvature  in  the  experimental  hysteresis  curves  (Figs.  4 and 
5)  at  zero  fields  is  an  obvious  feature  which  has  not  been  delineated  in  the 
present  calculations. 
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CAPTIONS 


1.  Structure  of  semicrystalline  polyvinyl idene  fluoride  showing  lamellar 
crystals  with  molecular  segments  normal  to  the  lamellae. 

2.  Assumed  potential  energy  of  a single  molecular  segment  rotating  about 
its  chain  axis  as  a function  of  angle  between  its  dipole  moment  and 

a direction  fixed  in  the  crystal.  The  fraction  of  all  molecules  in  a 
crystal  havinq  orientation  i is  given  by  f^ . 

3.  Calculated  polarization  hysteresis  for  the  6-site  model  assuming  lOnm 
long  molecular  segments  rotating  about  their  long  axis. 

4.  Experimental  polarization  hysteresis  for  a sample  cycled  at  increasingly 
higher  fields.  The  maximum  residual  polarization  (at  zero  field'  is 
about  lOuC/cm-. 

5.  Experimental  (Reference  3)  and  calculated  infra  red  transmission 
hysteresis  for  a vibration  polarized  along  the  CF-  dipole  in  the  s 
crystal  phase. 

6.  X-ray  pole  figure  data  from  Reference  3 (error  bars’*  and  calculated 
results  from  the  6-site  ferroelectric  model  (dasoed  line)  and  an 
equi librium  distribution  (solid  line). 
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Figure  2.  Assumed  potential  energy  of  a single  molecular  segment  rotating  about 
its  chain  axis  as  a function  of  angle  between  its  dipole  moment  and  the 
electric  field  direction.  The  fraction  of  uli  molecules  in  a crystal 
having  orientation  i is  given  by  f:. 
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Figure  6.  X-ray  pole  figure  data  from  Reference  5 (error  bars)  and  calculated 
results  from  the  6-site  ferroelectric  model  (dashed  lines)  and  an 
equilibrium  distribution  (solid  line). 
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